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Abstract 

We examine the possible consistency of the quasi-steady state model 
with the newly discovered SNe la. The model assumes the existence 
of metallic dust ejected from the SNe explosions, which extinguishes 
light travelling over long distances. We find that the model shows a 
reasonable fit to the data, which improves if one takes account of the 
weak gravitational lensing effect of the SNe which have been observed 
on the brighter side. 
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1. INTRODUCTION 

The high redshift supernovae (SNe) la explosions look fainter than they are 
expected in the Einstein-deSitter model, which used to be the favoured model 
before these observations were made a few years ago. This observed faintness 
is generally explained by invoking some hypothetical source with negative 
pressure often known as 'dark energy', the simplest and the most favoured 
candidate being a positive cosmological constant A. This happens because 
the metric distance of an object out to any redshift can be increased by incor- 
porating a 'fluid' with negative pressure in Einstein's equations. However, 
a constant A, is plagued with the so called cosmological constant problem: 
why don't we see the large vacuum energy density pv = A/SttG ~ 10^^ GeV^, 
required to drive inflation in the primordial epochs of the universe which is 
RilO^'^^ times larger than the value required by the SNe observations? It is 
more natural to believe that A dropped to zero, after the inflation was over, 
rather than try to explain its relic at such a small but highly fine-tuned value. 
Variable A or 'quintessence' models also fail to solve this problem without 
proper fine tuning, apart from the fact that they have a somewhat ad-hoc 
nature and certainly do not share the elegance of the overall structure of 
general relativity. 

An alternative way to explain the faintness of the high redshift SNe la 
is to consider the absorption of light by metallic dust ejected from the SNe 
explosions— an issue which, in the standard approach, is generally ignored 
while discussing m-z relation for SNe la. However, there is at least one theory, 
the quasi-steady state cosmology (QSSC), which considers this effect as a 
main ingredient of the theory. The QSSC is a Machian theory proposed by 
Hoyle, Burbidge and Narlikar (1993, 2000) as an alternative to the Standard 
Big Bang Cosmology (SBBC). This cosmology does not have any cosmic 
epoch when the universe was hot and is free from the initial singularity of 
the SBBC. It has been shown earlier (Banerjee, et al 2000; Vishwakarma 
2002; Narlikar, et al 2002) that by taking into account the absorption of 
SNe light by the intergalactic metallic dust, the QSSC explains successfully 
the SNe la data from Perlmutter et al (1999), together with SN 1997ff, the 
highest redshift SN observed so far (Riess, et al 2001). Extending that work 
further, we examine in this paper how well (or badly !) the QSSC fits the 
new data discovered with the Hubble Space Telescope (HST) which include 
7 highest redshift SNe la known, all at z > 1.25 (Riess, et al 2004). 

Being based on a Machian theory of gravity, the magnitudes and signs 
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of the creation field energy and A are determined by the large scale struc- 
ture of the universe (Hoyle, et al 1995). The role of A (which is negative in 
this theory) in the dynamics of the model is to energize the creation-field by 
controlling the expansion of the universe. The repulsive effects (akin to that 
from a positive A in the SBBC) are generated by the creation field which 
has a negative energy density. Therefore the model has cycles of expansion 
and contraction (regulated by the creation- and the negative A-fields respec- 
tively) of comparatively shorter period (around 50 Gyr) superposed on a long 
term (around 1000 Gyr) steady state-like expansion. Creation of matter is 
also periodic, being confined to pockets of strong gravitational fields around 
compact massive objects. These creation centers are 'turned on' close to the 
minimum scale size in a typical cycle and are gradually 'turned off' during 
expansion to maximum scale size. 

The theory seems to meet all the available observational constraints at 
the present time. According to the QSSC, the cosmic microwave background 
(CMB) is the relic starlight left by the stars of the previous cycles which has 
been thermalized by the metallic whisker dust emitted by the supernovae. It 
is very interesting to note that the energy available from this process is just 
right to give a radiation background of 2.7 K at the present epoch (Hoyle, et 
al 1994). SBBC, on the other hand, does not predict the present temperature 
of the CMB. The theory also explains the observed anisotropy of CMB, in- 
cluding the peaks at / ~ 200 and / ~ 600 which are related, in this cosmology, 
to the clusters and groups of clusters (Narlikar, et al 2003). The theory does 
not face the cosmological constant problem mentioned earlier. In fact, the A 
in the QSSC does not represent the energy density of the quantum fields, as 
this model does not experience the energy scales of quantum gravity except 
within the local centres of creation. 

It may also be noted that, unhke SBBC, the QSSC allows the dark matter 
to be baryonic. It may be recalled that the SBBC predicts the existence of 
non-baryonic, though as yet undetected, particles to solve the problems of 
structure formation and of the missing mass in bound gravitational systems 
such as galaxies and clusters of galaxies. Although there has been a steady 
evolution of views on whether dark matter is predominantly cold or hot, 
there is no satisfactory observational evidence for the postulated particles 
from laboratory physics. The predicted density distribution of dark halos 
which result from N-body simulations (Navarro, et al 1996), appears to be 
inconsistent with observations of spiral galaxies (de Blok, et al 2001) or with 
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strong lensing in clusters of galaxies (Treu, et al 2003). It is therefore fair 
to say that this scheme has still to demonstrate its viability. However, in 
the framework of the QSSC, the dark matter need not be necessarily non- 
baryonic. It can be in the form of baryonic matter being the relic of very old 
stars of the previous cycles. 

For the sake of completeness and for the ready reference, we describe 
briefly the mathematical formulation of QSSC in the Appendix. The details 
of this cosmology can be found in the papers of Hoyle, Burbidge and Narlikar 
mentioned above and in the paper by Sachs et al (1996). 

2. EXTINCTION BY METALLIC DUST 

Chitre and Narhkar (1976) were the first to discuss the role of inter galactic 
dust in the m-z relation, which was however largely ignored at that time. 
It is, however, generally accepted now that the metallic vapours are ejected 
from the SNe explosions which are subsequently pushed out of the galaxy 
through pressure of shock waves (Hoyle & Wickramasinghe 1988; Narlikar, et 
al 1997). Experiments have shown that metaUic vapours on cooling, condense 
into elongated whiskers of 0.5 — 1 mm length and falO"^ cm cross-sectional 
radius (Donn & Sears 1963; Nabarro Sz Jackson 1958). Indeed this type 
of dust extinguishes radiation travelling over long distances (Aguire 1999; 
Banerjee, et al 2000; Narhkar, et al 2002; Vishwakarma 2002; 2003). The 
density of the dust can be estimated along the lines of Hoyle, et al (2000). 
If the metallic whisker production is taken as 0.1 Mq per SN and if the SN 
production rate is taken as 1 per 30 years per galaxy, the total production 
per galaxy (of spatial density ~ 1 per 10^^ cm^) in 10^*^ years is ~ 2/3 x 10^^ 
g. The expected whisker density, hence, becomes 2/3 x 10^^ x 10"''^ ~ 10~^^ 
g cm"^. We shall see later that this value is in striking agreement with the 
best-fitting value coming from the SNe la data. 

Earlier work (Banerjee, et al 2000; Vishwakarma 2002; Narlikar, et al 
2002) has shown that the extinction due to dust adds an extra magnitude 
Amlz) to the apparent magnitude m{z) of a supernova of redshift z, where 

Am{z) = 1.0857 x k p^o (1 + z f-^rT^. (1) 

Jo ^[^) 

Here k is the mass absorption coefficient which is effectively constant over a 
wide range of wavelengths and is of the order 10^ cm^ g~^ (Wickramasinghe 
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& Wallis 1996); and pgo is the whisker grain density at the present epoch: 
Pg = Pgo Sq. The net apparent magnitude is then given by 

m''''\z) ^m{z) + Am{z). (2) 

The usual apparent magnitude m{z) arising from the cosmological evolution 
is given by 

m{z)^5log[Ho d^{z)]+M, (3) 
with the luminosity distance c^l given by 

for the k = case of the RW metric. The constant Ai appearing in equation 
(3) is given by = M — 5 logi/o + constant, where M is the absolute 
magnitude of the SNe. The Hubble parameter H{z) appearing in equations 
(1) and (4) is given by (A.16). 

3. DATA FITTING 

We consider the data recently published by Riess et al (2004) which, in 
addition to having previously observed SNe, also include 16 newly discovered 
SNe la by the HST, 6 of them being among the 7 highest redshift SNe la 
known, all at redshift > 1.25. We particularly focus on their 'gold sample' 
of 157 SNe la which is claimed to have a 'high confidence' quality of the 
spectroscopic and photometric record for individual supernovae. We note 
that the data points of this sample are given in terms of distance modulus 
fio = m"'^* — M = 51og(iL + constant. However, the zero-point absolute 
magnitude or Hubble constant were set arbitrarily for this sample. Therefore, 
while fitting the data, we can compare the observed Ho with our predicted 
j^net gjygn ]yy equatlou (2) and compute from 
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(5) 



i=l L "^Mo,! 



the constant M. thus playing the role of the normahzation constant. The 

quantity o"^^ . is the uncertainty in the distance modulus po,i of the i-th SN. 
We consider the simplest QSSC model with k = case of the RW metric 
(A. 4). Thus there are only three independent free parameters to be estimated 
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Figure 1: The allowed regions by the 'gold sample' of SNe la data (with 157 points) 
are shown at the 95% (inner contour) and 99% (outer contour) confidence levels, 
by marginalizing over A4. (The parameters k, pgQ and Hq are measured in units of 
10^ cm^ g~^, 10~^^ g cm~'^ and 100 km s~^ Mpc~^ respectively.) For an average 
values K = 5 X 10^ cm^ g~^ and Hq = 70 km s~^ Mpc~^, the estimated range of 
the whisker grain density from the data is 3.5 x 10~^^ g cm~'^ < pgo < 8 x 10~^^ 
g cm~^ at 99% confidence level. 

from the data: A^, kp^qIIq^ and f^Ao (only the last one comes from the field 
equations, see Appendix). The parameter kp^qIIq^, which is dimensionless, 
is of the order of unity if one considers k of the order 10^ cm^ g~^, pgo of the 
order 10"^"^ g cm~^ and Hq ~ 70 km s~^ Mpc~^. However, we have kept it 
as a free parameter to be estimated from the data. 

By varying the free parameters of the model, we find that decreases 
as f^Ao increases. Thus for the theoretically allowed region f^Ao < 0, the 
best-fitting is > 201.41 [at 155 degrees of freedom (dof), i.e., x^/dof = 
1.3]. For example, for the models 



= -0.1: xVdof = 205.85/155 
= -0.2: xVdof = 210.36/155 



1.33; 
1.36; 
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= -0.3: xVdof = 214.91/155 = 1.39. 

In order to compare these results with the SBBC, we note that for a constant 
A, the best-fitting fiat model (Omo + ^ao = 1) and the global best-fitting 
model (without any such constraint) are obtained as 

Qao = 1 - n^o = 0.69, with xVdof = 177.07/155 = 1.14; 
Qao = 0.98, Q^o = 0.46, with xVdof = 175.04/154 = 1.14. 

The so called ^concordance'' model Qmo — ^ — ^ao — 0.27 gives — 178.17, 
which is slightly higher than the best-fitting value of for the flat model 
mentioned above. In fact, there is an almost fiat valley around f^mo = 0.3 
on the Q^-ao-Ai-x"^ surface where x"^ does not vary significantly and hover 
around 177-178. For example, the models f2mo(= 1 — ^ao) = 0.28 and 0.3, 
respectively, give x^ — 177.67 and 177.13. 

Though the fit in the QSSC is certainly not as good as in the SBBC, it is 
by no means rejectable. Moreover the estimated value of kp^qH^^ is indeed 
of the order of unity: the best-fitting values of kp^qHq^ for the cases VLao = 
-0.1, -0.2 and -0.3 are respectively 4.19, 4.75 and 5.30, as expected from the 
theory. Here, the parameters k, pgo and Hq are measured in units of 10^ cm^ 
g"^, 10~^^ g cm~^ and 100 km s~^ Mpc~^ respectively. In Figure 1, we have 
shown the allowed regions in the parameter space ^aq — kp^qHq^ at 95% 
and 99% confidence levels, by marginalizing over the parameter M.. Figure 
2 compares the best-fitting theoretical models with the actual data points. 

4. EFFECTS OF WEAK LENSING 

Weak gravitational lensing is an unavoidable systematic uncertainty in the 
use of SNe la as standard candles. As the universe is inhomogeneous in mat- 
ter distribution, the SNe fluxes are magnifled by foreground galaxy excess 
and demagnified by foreground galaxy deficit, compared to a smooth matter 
distribution. Recently Williams & Song (2004) have reported such a correla- 
tion between the magnitudes of 55 SNe from the sample of Tonry et al (2003) 
and foreground galaxy overdcnsity. They have found the difference between 
the most magnifled and the most demagnifled SNe as about 0.3-0.4 mag. 
Wang (2004) has claimed further evidence of gravitational magniflcation of 
three brightest SNe from the Riess et al (2004) sample: 

SN1997as {z = 0.508, /x« = 41.64): magnified by 2.10 ± 0.68; 
SN2000eg {z = 0.540, Po = 41.96): magnified by 1.80 ± 0.70; 
SN1997as {z = 0.886, po = 42.91): magnified by 2.42 ± 1.98. 
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Figure 2: Best-fitting flat models in QSSC and SBBC are compared with the 'gold 
sample' of SNe la data from Riess et al (2004). The QSSC model corresponds to 
f^AO = -0.2. 

These high magnification factors ~ 2 from weak lensing are though somewhat 
surprising, as also mentioned by Menard & Dalai (2004) who claim not to 
find any significant correlation between SN magnification and foreground 
galaxy overdensity. However, even if we assume a mild magnification of 
the above-mentioned SNe just by an average 0.5 mag, this improves the fit 
considerably: 



^AO = 


-0.1: 


-^improved /'^^■^ 


198.55/155 = 


1.28; 


^AO = 


-0.2: 


-^improved /'-^^■^ 


203.17/155 = 


1.31; 


^AO = 


-0.3: 


Ximproved/'^^-^ 


207.82/155 = 


1.34. 



Though the weak lensing effects are estimated to be small for SNe at z < 1, 
they are non-negligible for higher redshift SNe. As more SNe are discovered 
at higher redshifts, it becomes increasingly important to minimize the effect 
of weak lensing, say, by considering the flux-averaging (Wang, 2000). 

5. CONCLUDING REMARKS 
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The redshift magnitude test has had a chequered history. During the 1960s 
and the 1970s it was used to draw very categorical conclusions. The deceler- 
ation parameter qq was then claimed to lie between and 1 and thus it was 
claimed that the universe is decelerating. Gunn and Oke (1975), however, 
pointed out that there remained observational errors to be allowed for, that 
vitiated that conclusion. For example, aperture correction, luminosity evo- 
lution, etc. were to be allowed for. It was then realized that the test is not 
as conclusive in selecting a cosmological model as it earlier appeared to be. 

Today's situation, we feel, is hardly different. There has been considerable 
progress in our understanding of the physics of supernovae, yet it is hard to 
imagine that the peak luminosity of Type la supernovae remains a standard 
candle over a redshift exceeding 1. Evolution has long been assumed in 
various other cosmological tests, like the counts of galaxies and radio sources, 
the variation of angular size with redshift, etc. The assumption of a non- 
evolving standard candle for supernovae therefore needs to be more critically 
examined than has been hitherto. 

The possible role of gravitational lensing in amplifying the supernova 
luminosity at high redshifts has been discussed by several authors and we 
have applied those ideas here to illustrate the difference it can make to any 
conclusion drawn from the data. Additionally, the role of intergalactic dust 
still remains to be appreciated fully and we have demonstrated here the 
possible difference it can make to the viability of a model. 

Contrary to the widespread belief that these caveats do not matter or 
have already been allowed for, we retain a healthy skepticism of this test as 
contributing to a 'precise' determination of cosmological parameters. For this 
reason we arc satisfied with the level of 'goodness of fit' obtained here for the 
QSSC. The fit could no doubt be improved by tinkering with the parameters; 
but given the observational uncertainties, we do not feel it worthwhile to 
undertake that exercise. 
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APPENDIX 

The field equations of QSSC, which arise from a Machian theory of gravity, 
are more general than the Einstein field equations: 

RiJ _ g^J _ A = _87rG [T^,tter + T^Ltion] • (A-l) 

Here the speed of light is taken as unity. The first term on the right is the 
usual energy momentum tensor of matter 

T^.tter^ip + pyu^+pg'^, (A.2) 

whereas the second term denotes the contribution from a trace-free zero rest 
mass scalar field c of negative energy and stresses with gradient q = dc/dx^: 

^cLtion = -/ (cV + ^C^Q g'^^ , (A.3) 

with a positive coupling constant /. In the case of the homogeneous isotropic 
spacetime described by the RW metric 

ds' = -dt' + S\t) I + r\de' + sin^ e d<P') I , (A.4) 

the field equations (A.l— A.3) lead to the following two equations for the 
'dust' universe {p — 0): 

k A 8nG. , 3^.2 

■g^+g^ = ^ + ^iPr^ + Pc), Pc = -^fc] (A.5) 

2| + |^ + ;^ = A-87rGpe, P.^-\fc?. (A.6) 

By assuming that the present epoch is represented by the non-creative mode 
of the model, i. e., T^^^^ter -j = T'^e^tion -j = 0, giving 

8nG A 8nG B ^ ^ , , 

-^Pm = -^Pc = A,B = constants, (A.7) 

equation (A.5) can be solved to give the scale factor in the form 

S = S[l + 7] cos ^lj{t)]. (A.8) 
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The parameter rj lies in the range < < 1 and the function ■0 is given by 

ip^ ^ --(l+r7cos^)"^{6 + 4r7cosV' + ^?^(l + cos^^)}. (A. 9) 

3 

Obviously A is negative, as has been mentioned earlier. It is also obvious 
from equation (A. 8) that S never becomes zero and oscillates between 

S{1 -rj) = Srnin <S<S^^ = S{1 + rj). (A.IO) 

The constant S, appearing in (A. 8) and (A.IO), is given by 

A^2kS -^AS^{l + r]^), (A.ll) 

B = kS'il - rf) - \kS\1 - r,'){3 + v'), (A.12) 

which can be obtained from equations (A. 5) and (A. 6) using (A. 7) and (A. 8) 
therein. Equations (A.ll) and (A.12) can be recast in the following forms in 
terms of the different energy components computed at the present epoch: 

J^mo = 2f]fco(l + z)-^ - 4f^Ao(l + z)-%l + T}^), (A.13) 

^^co = -^^ifco(l + z)-^(l - ri") - 4Qao(1 + z)-\l - ri'){?> + rf), (A.14) 
where, as usual. 



o - o - ^ r, - A o - ^ A 1 



In terms of these dimensionless parameters, equations (A. 5) and (A. 6), by 
the use of (A. 7), reduce to 

H{z) = i/o[^^Ao - ^^feo(l + zf + 0^o(l + zf + ao(l + -2)1'/', (A.16) 



2q{z) 



[^^mo(l + z)' - 20ao + 20,o(l + z)% (A.17) 



H{z)_ 

It should be noted that not all the parameters, introduced above, are inde- 
pendent. For example, equation (A.16) suggests that 

^AO ~ ^fcO + ^mO + ^cO = 1- (A. 18) 
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This equation also suggests that at the maximum redshift 2;max (say, in the 
present cycle), one has the following identity: 

QaO - Ofeo(l + Z^^f + 0„,o(l + Z^^f + Oeo(l + -2max)^ = 0. (A.19) 

Also equation (A. 10) suggests that z^ax., z and r) are related by 

l + z = ^^{l-r)){l + w). (A.20) 

Thus out the 7 parameters Qao, f^mO; ^cO; ^ko, z, rj and -Zmax, only 3 parame- 
ters, say, QjfcO) ^Ao and ^r^ax are independent. For the case A; = 0, which gives 
the simplest one of the QSSC models, only two independent parameters Qaq 
and Zmax are left out. We further consider z^-aax = 5, as has been done in the 
earlier papers on the QSSC, which leaves only one free parameter, say, f^Ao 
coming from the field equations. 
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